Experimental protection and revival of quantum correlation in open solid systems 
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Quantum correlation quantified by quantum discord has been demonstrated experimentally as 
important physical resources in quantum computation and communication for some cases even 
without the presence of entanglement. However, since the interaction between the quantum system 
and the noisy environment is inevitable, it is essential to protect quantum correlation from lost 
in the environment and to characterize its dynamical behavior in the real open systems. Here 
we showed experimentally in the solid-state P:Si system the existence of a stable interval for the 
quantum correlation in the beginning until a critical time t c ~ 166 ns of the transition from classical 
to quantum decoherence. To protect the quantum correlation, we achieved the extension of the 
critical time by 50 times to 8 ps by applying a two-flip dynamical decoupling (DD) pulse sequence. 
Moreover, we observed the phenomenon of the revival of quantum correlation, as well as classical 
correlation. The experimental observation of a non-decay interval for quantum correlation and the 
great extension of it in an important solid-state system with genuine noise makes the use quantum 
discord as physical resources in quantum information processing more practicable. 
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Quantum correlations that can be shared only among 
quantum systems play essential roles in fundamental 
physics and future technologies p]. Although quantum 
entanglement has been widely investigated, it cannot de- 
scribe the nonclassicality of the correlations in separable 
states [2J, while quantum discord can [3J[3]. A great deal 
of theoretical investigations [ST-|13| on quantum discord 
have been carried out recently to show the properties 
and importance of quantum discord. Moreover, quan- 
tum discord has been demonstrated experimentally as 
important resources in deterministic quantum computa- 
tion with one pure qubit (DQC1) |14) and remote state 
preparation even without the presence of entanglement 
|15j . Furthermore, an operational method to use quan- 
tum discord (without the presence of entanglement) as 
a physical resource has been demonstrated in Ref. |16| . 
Since all quantum systems inevitably interact with the 
environment, the evolution of the quantum discord in a 
noisy environment is certainly of great interests. It is 
well-known that in a dissipative environment entangle- 
ment decays to zero in a finite time, which is called entan- 
glement sudden death (ESD) [TB]. In contrast, quantum 
discord has been claimed to have no sudden death [T§] . 
but has peculiar sudden change in its decay rates [20] . It 
has also been shown that quantum discord is robust in its 
initial period of decoherence but suffers a sudden change 
phenomenon |21) . Such a behavior has been recently ex- 
perimentally investigated in optics [55] and liquid NMR 
systems |23j . but the dynamics of quantum discord in 
solid-state systems with a real noisy environment remains 
elusive. 

Meanwhile, it is essential to preserve the quantum cor- 
relation in a fragile quantum system, thus it is vital to 
overcome decoherence effect induced by the environment. 



One of the strategies is the dynamical decoupling (DD) 
technique [2H [25] , which uses stroboscopic spin flips to 
reduce the average coupling to the environment to ef- 
fectively zero. DD is a particularly promising strategy 
for combating decoherence [21H |2"T] , since it can be nat- 
urally integrated with other desired functionalities, such 
as quantum gates. 

Herein, we investigate the dynamics of both the clas- 
sical and quantum correlations in solids. We find that 
in the initial interval, the quantum correlation is stable 
while the classical correlation decays. After this time 
interval, the quantum correlation decays, while the clas- 
sical correlation remains. Furthermore, we manage to 
extend the time point of the sudden transition from the 
classical decoherence to the quantum decoherence is de- 
layed from about t c sa 166 ns to about 8 [is by DD. This 
time point can be further extended to about 12 /is by 
a multiple two-flip DD pulse sequence. Additionally, we 
have observed the revival of both the classical and the 
quantum correlations, which may have some impact on 
quantum information processing. 

For a biparticle state pab, the total correlation is quan- 
tified by mutual information I {pab) — S{pa) + S{ps) — 
S{pab), where S{p) = — Tr[pLog 2 p] is the von Neumann 
entropy, pa and ps are the reduced-density matrices of 
Pab- When performing a positive-operator- valued mea- 
sure {11^} on particle B, the maximal accessible infor- 
mation about particle A gives the classical correlation 
C b {pab) = max {Tl B } [S{pA) - J2kPkS{p A )}, where p k = 
Tr[(l®nf) PAB (l®ITf)t], p\ = Tr s [(l®nf) PAB (l® 
n^)i]/|?fc. Hence, the quantum part of the total correla- 
tion is defined as quantum discord Db{pab) = I{pab) — 
Cb{pab)- Note that quantum discord is asymmetric, 
i.e. generally Db{pab) 7^ Da{pab)- Herein, for sim- 
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plicity we consider an initial states of the form of the 
Bell-diagonal states, 

PAB = \{i+Y,avt°f), a) 

i=l 

where crf^ B '(i — 1,2,3) are the Pauli matrices. For the 
density matrices of the Bell-diagonal states, Ref. [25] and 
[29] provided the analytical expression and the geometric 
picture of the quantum discord, respectively. 

In our experiment, we have chosen phosphorous donors 
in silicon (P:Si) [27J [3T] [33] material with P concentra- 
tion about 1 x 10 16 cm~ 3 for the study. Silicon is a par- 
ticularly attractive material for hosting spin qubit for its 
low-orbit coupling and low natural abundance of nuclear- 
spin-bearing isotope |30j . The experiment was performed 
at temperature 8 K. P:Si consists of an electron spin, 
S = 1/2 (<? = 1.9987), coupled isotropically to the nu- 
clear spin, / = 1/2, of 31 P, with a hyperfme coupling 
constant A = 117 MHz. The energy diagram of this 
system is plotted in FigJTJi, where the four-level system 
can be manipulated by resonant the microwave (MW1, 
MW2) and radio- frequency (RF1, RF2) radiation. The 
dynamics of an open quantum system pab coupling to 
the environment in solids can be described by longitudi- 
nal and transverse relaxations of time constants T± and 
T 2 , respectively. For the electron spin, the transverse 
relaxation time is T 2e — 120 fis, the longitudinal (i.e., 
electron population) relaxation time is T± e = 5.6 ms and 
the dephasing time T 2e ~ 0.2 /is. For the nuclear spin, 
the dephasing time is determined as T 2n = 24 ps using 
the nuclear spin free induction decay (FID) experiment. 

The longitudinal relaxation time T\ of both the elec- 
tron and the nuclear spins are much larger than the trans- 
verse relaxation time T 2 and so can be neglected in our 
experiments. The decay of the off-diagonal elements de- 
pends on the time scale of dephasing time T 2 of the elec- 
tron and the nuclear spin. Since the electron dephasing 
time is almost two orders of magnitude smaller than the 
nuclear dephasing time, the decay of the secondary diag- 
onal elements is dominated by T 2e (see Supplementary 
Information for details). For the phosphorous donors 
in silicon with natural abundance of 29 Si, the hyperfme 
fields of the 29 Si nuclei cause random static shifts of the 
individual electron-spin resonant frequencies which sat- 
isfy a Gaussian distribution [33]. Averaging over the 
Gaussian-distributed resonant frequencies, the dephas- 
ing of the electron spin can be derived with the form 
~ exp[— (i/T 2 * e ) 2 ]. As a result, the time evolution of the 
total system is given by 

PAB {t) = ^[l+ Cl (t)^CT^+c 2 (t)a^a^+c 3 (t)a^af}, 

(2) 

where the coefficients Ci(i) = ci(0) exp[— (t/T 2e ) 2 ], 
c 2 (t) = c 2 (0) exp[-(t/T| e ) 2 ] and c 3 (t) = c 3 (0) = c 3 . 
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FIG. 1. Energy level diagram, experimental pulse se- 
quence and reconstructed deviation density matrix. 

a, Energy level diagram for the P:Si system. There are four 
Zeeman product states which are labeled by states 1 — 4, 
respectively, f and 4- stand for the ±1/2 states of electron 
and nuclear spins. EPR and NMR transitions are indicated 
by two-way arrows, b, Diagram of the experimental pulse 
sequence, which includes three steps: initial state prepara- 
tion, relaxation delay and final state detection, c, d, The 
real and imaginary parts of the reconstructed deviation den- 
sity matrix Ap — p — \, in unit of e, respectively, where 
s = g/3 e B /8k B T = 7.35 x 1(T 3 (at temperature 8 K) is the 
ratio between the magnetic and thermal energies. 

The typical coefficients Cj initialized by electron spin 
resonance (using ~ 10GHz excitation at a temperature of 
8 K) are ~ 10~ 3 , therefore we focus on a class of reason- 
able and widely used states for which ci(0) = 0, |c3(0)| < 
|c 2 (0)| <C 1. By following Ref. [35], we obtain the analyt- 
ical expressions for the mutual information, the classical 
correlation and the quantum discord. Expanding I[p(t)] 
and C[p(t)] in the Taylor series of c 2 (£) and C3 and ne- 
glecting high-order terms, we obtain 

i[p(t)] = ^[4 + 4(t)}, (3) 



C[ P (t)] = ^2 c2 W' W 

where c(t) = max{c 2 (i), C3}. Hence, the quantum discord 
is calculated to be 

DW)] = V^^ (5) 

Here t c = — ln[c 3 /c 2 (0)]T 2e is obtained by setting 
ciitc) — C3. Consequently, in the noisy environment, the 
quantum discord is constant and the classical correlation 
decreases in the initial period t ^ t c , while for t > t c , 
the classical correlation does not change in time and only 
the quantum discord is reduced sharply. 

In our experiment, the duration of the microwave and 
the radio-frequency n pulses are 60 ns and 10 ps, which 
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FIG. 2. Dynamics of classical and quantum Correla- 
tions. The logarithmic-scale values of mutual information 
(darkgreen circle), quantum discord (red square) and classi- 
cal correlation ( black triangle) are numerical computed with 
their original definitions (and so do the other experimental 
data in the article) (see Supplementary Information for de- 
tails). The curves are the theoretical predication according 
to Eq. JjjJ. I and II stand for the classical decoherence and 
quantum decoherence regimes, respectively. Error bars are 
calculated (see Supplementary Information for details). 



are determined using the EPR and NMR Rabi nutation 
experiments, respectively. By following Ref. [3T], the nu- 
clear population relaxation time T\ n is estimated to be 
250 times of T\ e . So the waiting time between each ex- 
periment is set to 10 s. FigjlJ) shows the first pulse se- 
quence applied in our experiment. It consists of three 
steps: the initial state preparation, the relaxation delay 
and the final state detection. Starting from the thermal 
equilibrium state, a MW2 pulse is used to flip the electron 
spin with an angle 9\ while the nuclear spin I z = 1/2. 
This is followed by a waiting time of (ti = 1 fj, s 3> T^), 
so that the off-diagonal elements of the density matrix 
decay off. After that the tt/2 R- F2 (RF1) pulse and the 
ir MW2 pulses are applied to generate the nonzero sec- 
ondary clement. Since the pulses tt/2 RF1 and 7r/2 RF2 
are applied sequentially, after generating the off-diagonal 
elements P34 and P43 by RF2, they will decay during 
the period of the next pulse RF1 with T| n . To compen- 
sate for this effect, we add a 6*2 RF1 between the first 
pulse MW2 and the duration r 2 (r 2 = 200 ps > T|J, as 
shown in the FigfTJ) (see Supplementary Information for 
details). In our experiment, the actual values of 9\ and #2 
are chosen to be 6q w 0.70 7r and 62 ~ 0.28 ir. Then the 
initial state is prepared as a Bell-diagonal state, which 
satisfies the conditions ci(0) = 0, | C3 (0) | < |c2 (0) | <C 1. 
We present the result of the tomography of the initial 
state in Fig{TJ;,d, which represent the real and imaginary 
parts of the deviation density matrix in unit of e. Af- 
ter the preparation step, the quantum system is left to 
evolve under the noisy environment. Then the resultant 
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FIG. 3. Preservation of classical and quantum correla- 
tions by dynamics decoupling, a, Schematic illustration 
of the dynamics decoupling pulse sequence. After the initial 
state preparation, a set of two-flip DD sequence pulses for the 
electron spin are applied at r and 3r, and then the state to- 
mography is performed at 4r. b, The plot of the results using 
the same notation as in Fig. 2. The curves are drawn to follow 
the trend of the variation of correlations. I and II denote the 
classical and quantum decoherence regimes, respectively. 



quantum states are reconstructed by state tomography 



In our case, according to Eqs. Q and ([5]), the dynam- 
ics of the quantum and classical correlations depend on 
the decay of secondary diagonal elements. The decay 
of P23 an d P14 is shown in Supplementary Information 
with the fitted decay function of ~ exp[— (i/^decay) 2 ], 
where Td ccay = 175 ns. The values of mutual information 
(darkgreen circle), the quantum discord (red quadrate) 
and the classical correlation (black triangle) calculated 
from the state tomography results (see Supplementary 
information for details) are plotted in Figj2] The theo- 
retical predications plotted in Figj2] can successfully de- 
scribe the experimental data and clearly show that the 
sudden transition from the classical to the quantum de- 
coherence regime occurs at about 166 ns. In the initial 
period (0 < t < 166 ns), the quantum discord remains 
constant but the classical correlation decreases, while for 
t > 166 ns, the classical correlation does not change but 
the quantum correlation decreases dramatically. 

Dynamics decoupling is an effective method to pre- 
serve the coherence of a quantum system, and it has 
been tested by experiments in both single qubit |26] and 
two-qubit [27] systems. Herein we use the two-flip DD se- 
quence pulses for the electron spin to prolong both classi- 
cal and quantum correlations. Fig|3^, shows the sequence 
pulses used in the dynamics decoupling experiment. Af- 
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FIG. 4. Revival of classical and quantum correlations 
by dynamics decoupling pluses, a, Schematic illustration 
of the revival experiment pulse sequences. After the initial 
state preparation, multiple two-flip DD pulses (r^ = 1 /is) for 
the electron spin are applied, b, The plot of the results in 
logarithmic scale using the same notation as in Fig. 2. I and 
II stand for the classical decoherence and quantum decoher- 
ence regimes, respectively. Ill and IV stand for the quantum 
revival and classical revival regimes, respectively. 

ter preparing the initial state, the two-flip DD sequence 
pulses for the electron spin are applied, and the finally 
state tomography is performed at At. FigjSjD shows that 
the decay of classical and quantum correlations become 
much slower, and the period before the sudden transition 
of correlations is prolonged via DD by about 50 times. 

After the dephasing, the classical and quantum corre- 
lations decay to almost zero, but do not disappear com- 
pletely. When a two-flip DD sequence pulse for the elec- 
tron spin is applied after a period T4, as the Fig|4£i shows 
(T4 = 1 /is), along with another same two- flip DD se- 
quence pulse for the electron spin at time 3t4, there will 
be a revival of the classical and the quantum correla- 
tions at about the time 4t4. This process could be re- 
peated several times, and it is observed that both the 
classical and the quantum correlations vanish and then 
revive. Fig|4jD shows the experimental observation, it is 
quite clear that the quantum discord (red square) decays 
from a constant value to nearly zero and then revive to 
a stable value at different stages by applying a proper 
two-flip DD pulse sequence, while the classical correla- 
tion (black triangle) decreases dramatically from an ini- 
tial value to a constant value and has a similar revival. I 
and II stand for the classical decoherence and quantum 
decoherence regimes, respectively. The sudden transi- 
tion from a quantum revival regimes (labeled by III ) to 



a classical revival regimes (labeled by IV) is obtained in 
Fig. |1Jd. Moreover, we can see that the time interval be- 
fore the sudden transition of the correlation is prolonged 
to about 12 /xs by the sequence of the three two-flip DD 
pulses, longer than the time by applying a single two-flip 
DD showed in Fig|3|D. 

We are aware of that there are other measures for 
quantum correlations |35j . among which is the geomet- 
ric measure of the quantum discord [36 , which differs 
in general from the conventional definition of quantum 
discord ^ . This geometric measure of quantum dis- 
cord has been shown [15] to be related to the fidelity of 
quantum remote state preparation. For the two-qubit 
system studied here, the dynamical behavior of the ge- 
ometric measure of quantum discord turns out to be in 
coincidence with the conventional quantum discord (see 
Supplementary Information for details). 

We have studied experimentally the dynamical behav- 
ior of the quantum discord (quantum correlation) , as well 
as the classical correlation of P:Si system, a solid-state 
system with genuine rather than simulated noisy envi- 
ronment. The results show clearly a transition from clas- 
sical decoherence to quantum decoherence at a critical 
point t c ss 166 ns: in the initial interval < t < t c , the 
quantum discord does not change while the classical cor- 
relation decays, but after the critical point, the classical 
correlation stays constant, while the quantum discord de- 
cays. Furthermore, we have not only observed that the 
critical point t c can be extended via a simple DD ap- 
proach by about 50 times to 8.0 /is and further to about 
12 /xs with multiple two-flip DD pulses, but also clearly 
demonstrated the decay and revival of both the classical 
and the quantum correlations. In view of that quantum 
discord could be the reason for the power of quantum 
computation in some cases [SJ [H] , a resource for remote 
state preparation |15j . together with the demonstration 
of an operation method to use quantum discord as a phys- 
ical resource |16j . here the experimental demonstration 
of the existence of a non-decay region, the revival and 
prolonging of the quantum discord in a noisy solid-state 
system may have great potential applications in quantum 
information processing. 
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SUPPLEMENTARY INFORMATION FOR 
"EXPERIMENTAL PROTECTION AND 
REVIVAL OF QUANTUM CORRELATION IN 
OPEN SOLID SYSTEMS" 

1 STATE PREPARATION 

In order to prepare a density matrix in the form of the 
Bell-diagonal states which satisfy the conditions C\ (0) = 
0, |c3 (0) | < | C2 (0) | <C 1, several different microwave and 
radio-frequency pulses need to be designed elaborately. 
Starting from the thermal equilibrium state po = \ l4 X 4— 
ecr z ® 1 2X 2, where e = g(3 e B /8k B T = 7.35 x 1CT 3 (at 
temperature 8 K) is the ratio between the magnetic and 
thermal energies, a MW2 pulse is used to flip the electron 
spin with an angle 9\ while I z — 1/2. There is a waiting 
time (ti = 1 fj,s ^> T 2 * e ) followed to let the off-diagonal 
elements of the density matrix decay off. Then a 6* 2 RF1 
is applied, which is followed by t 2 = 200 p,s ^> T 2 * e . It 
is noted that, since both the electron and the nuclear 
spin longitudinal decoherence time Ti's are sufficiently 
larger than not only the transversal docoherence time 
T 2 , but also the experiment time scale, the longitudinal 
relaxation can be neglected, and so the diagonal elements 
of the density matrix remain invariant. 
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The next 7r/2 pulse RF2 is to equalize the diagonal 
elements P33 and P44, along with generating some other 



equivalent off-diagonal elements p 3 4 and P43. The coher- 
ence elements decay with a characteristic time T 2 * n . 
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The following ir/2 pulse RF1 is to equalize the diagonal 
elements pu and p 22 , along with generating some other 
equivalent off-diagonal elements pi 2 and p 2 i- Since the 



pulse length of the radio- frequency 717 '2 pulse is 5 /us, the 
decay of P34 and P43 in the period of pulse 7r/2 RF1 needs 
to be considered. 
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Where /(i§) is the proportion of ^34 and P43 decay 
during the 7r/2 pulse RF1. The last MW2 pulse is applied 



to transfer the off-diagonal elements P34, P43 to pu, pu, 
and the off-diagonal elements p\2, P21 to P23, p?,2- 



P4= 4 1 



/ cos 2 [^] 

-cos 2 [\\ cos[0 2 ]sin 2 {\\ 

cos[6» 2 ]sin 2 [\\ 

V-sin 2 [^]/(i f ) 



-sin 2 





-cos 2 



COS 





2 [0 



(9) 



m j 



A suitable choice of the angles 6\ and 2 satisfy- 
ing the conditions cos[0 2 ] = f(t^) and cos 2 [^] < 
cos[6* 2 ]sin 2 generates the final density matrix in the 



form of the Bell-diagonal state, which satisfies the condi- 
tions Ci(0) = 0, 1 c 3 (0) I < |c 2 (0)| <C 1. In our experiment, 
the actual values of 9\ and 9 2 are chosen to be 9\ « 0.70 ir 
and 62 ~ 0.28 n. So the theoretical density matrix after 
the five pulses becomes 



P5 = 7 1 



/ 0.206 -0.506\ 

-0.206 0.506 

0.506 -0.206 

\-0.506 0.206 j 



(10) 



STATE TOMOGRAPHY AND THE DECAY OF 
OFF-DIAGONAL ELEMENTS 



raphy is needed throughout the state evolution. The 



To identify the dynamics of classical and quantum cor- 
relations, besides the state preparation, the state tomog- 



7 



state tomography of two qubit in solid is used and in- 
troduced in detail in [3 2]. Here in our experiment, the 



density matrix of is measured after the state preparation 
to be 



PAB 



( 0.21(3) 0.01(1) + £0.02(1) 
0.01(1) - £0.02(1) -0.21(3) 

-0.01(1) + £0.02(3) 0.50(8) + £0.15(2) 

\-0.54(8) - £0.11(1) 0.02(1) - £0.01(1) 



-0.01(1) -£0.02(3) 
0.50(8) -£0.15(2) 

-0.17(3) 
0.01(2) +£0.02(2) 



-0.54(8) +£0.11(1)^ 
0.02(1) +£0.01(1) 
0.01(2) - £0.02(2) 
0.17(3) / 



(11) 
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FIG. 5. The decay of the real part of pn and p23- 
a, b, The decay of the real part of pi4 and P23, in unit 
of e, respectively. The solid line is the fitted function of 
~ exp[-(t/T dccay ) 2 ], with T dccay = 175 ns. 



The quoted errors come from the fitting errors for the 
nuclear nutation. 

According to Eq.(4) (main text) and Eq.(5) (main 
text), the discord and classical correlation of the states 
which we are studying just depend on the decay of /?23, 
pn and their relationship with the diagonal elements. 
The decay of p23 and pu are measured as Figj5^,b shows, 
in unit of e, respectively. The solid line is the fitted func- 
tion of ~ cxp[-(t/T decay ) 2 ], with T decay = 175 ns. 



tive for two-qubit states, it is enough to maximize over 
all the following projective measurements {1<E> | Ofc >< 
Q k |,fc=||,_L}, where | }l ) = cos (9 | 0) +e^sin0|l) and 
|0_i_) = e - ^ sin#|0) — cos#|l) presents an arbitrary ba- 
sis of B formed by two orthogonal states on the Bloch 
sphere, with < 9 < | and < <j> < 2tt. The classical 
correlation and quantum discord are calculated numer- 
ically to be 1.8(6) x 10~ 4 and 2(1) x 10~ 5 . The errors 
of the correlations are obtained with the maximum dif- 
ference of correlations between density matrices p' AB and 
PAB, where each element of p' AB falls near that of pab 
within the element errors. 



4 DYNAMICS OF GEOMETRIC MEASURE OF 
QUANRUM DISCORD 

There are many ways to quantify and verify quan- 
tum correlations [35] > such as quantum deficit, distillable 
common randomness, measurement-induced disturbance, 
symmetric discord, relative entropy of discord and disso- 
nance and geometric measures. Recently, it has been 
shown that the geometric measure of quantum discord is 
related to the fidelity of quantum remote state prepara- 
tion, which provides its operational meaning |15j . Herein 
we focus the dynamics of geometric measure of quan- 
tum discord which first introduced by Dakic, Vedral, and 
Brukner [36] . This measure is significant in capturing 
quantum correlations from a geometric perspective and 
can be evaluated explicitly and leads to an explicit for- 
mula for any two-qubit state. It is defined as the normal- 
ized trace distance to the set of classical states 15, 36 



3 CALCULATION OF CORRELATIONS 

Here we take the reconstructed density matrix of 
Eq.( |ll[ ) as an example, i.e. the one shown in Fig. lc,d 
(main text). The mutual information of pab is calcu- 
lated with I(pab) = S(p A ) + S(p B ) - S(pab) to be 
2.0(6) x 10~ 4 . The classical correlation and quantum 
discord are obtained with Cb(pab) = max {n^} [S(pa) ~ 

J2kPk s (PA)} and d b(pab) = I(pab) - C b (pab) by op- 
timizing over all 1-qubit measurements. Since it has been 
proven that the optimal measurement is always projec- 



V 2 (pab) = 2rmn xe n \\p A B-x\ 



= 2min xe o Tr( / 04 £r 
(12) 



where £Iq denotes the set of zero-discord states. 

For evolutive Bell-diagonal states Eq.(2) in the main 
text, the geometric measure of quantum discord can be 
expressed as 



V 2 W)] - j (^(*)+ci(*)+^(t)-niax{c?(t),^(*),ci(t)}), 

(13) 
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FIG. 6. Dynamics of quantum discord and geometric 
measure of quantum discord, a, Quantum discord (red 
square) and geometric measure of quantum discord (dark- 
green triangle) are numerical computed with their original 
definitions. The curves are the theoretical predication accord- 
ing to Eq.(5) (main text) and Eq.( |14[ ). I and II stand for the 
constant and decoherence regimes, respectively, b, Preserva- 
tion of quantum discord and geometric measure of quantum 
discord by dynamics decoupling. The curves are drawn to 
follow the trend of the variation of correlations. 



U <4(t) ift>t c 



(14) 



same to the main text t r 



In the case of ci(0) = 0, |c 3 (0)| < |c 2 (0)| <C 1, it is calcu- 



Where t c 

y-ln[c 3 /c 2 (0)]r 2 * e . Comparing Eq.([l4| with Eq.(5) 
(main text), we see that the difference between geometric 
measure of quantum discord and quantum discord is just 
a constant coefficient l/ln2. Consequently, in the noisy 
environment, geometric measure of quantum discord will 
be constant in the initial period t ^ t c , while after the 
period t c , it will be reduced sharply. 

Figjf^i shows the dynamics of quantum discord and 
geometric measure of quantum discord. The values of 
Quantum discord (red square) and geometric measure 
of quantum discord (darkgreen triangle) are numerical 
computed with their original definitions. There is sud- 
den transition of quantum discord from constant to de- 
coherence regimes occurs at about 166 ns as well as geo- 
metric measure of quantum discord. In the initial period 
(t < 166 ns), both the quantum discord and geometric 
measure of quantum discord remain constant, while for 
t > 166 ns, they decrease dramatically. Figj6] shows that 
the decay of both quantum discord and geometric mea- 
sure of quantum discord become slower, and the period 
before the sudden transition of correlations is prolonged 
by DD from about 166 ns to 8 [is. 



